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A special relativistic perturbation to non-relativistic quan- 
tum mechanics is shown to lead to the special relativistic pre- 
diction for the rate of precession for quantum states in the 
Coulomb potential. This behavior is illustrated using SO (4) 
coherent states as examples. These states are localized on 
Kepler ellipses and precess in the presence of a relativistic 
perturbation. 



I. INTRODUCTION 

Classical objects which are bound by forces which vary 
by the inverse square of distance have elliptical orbits 
whose orientation does not change with time. It is well 
known that the inclusion of relativistic effects causes 
these elliptical orbits to precess. In the Coulomb poten- 
tial, where interactions are electromagnetic, the rate of 
this precession is calculated using special relativity. This 
was first accomplished by Sommerfeld (lj who, for the 
purpose of calculating the fine structure of the hydrogen 
atomic spectrum in the old quantum theory, calculated 
the classical rate of precession due to a special relativis- 
tic treatment of the kinetic energy of a "spinless" elec- 
tron the Coulomb potential. (Sommerfeld's calculation 
ignored electron spin in part because it was not described 
until 1925 by Uhlenbeck and Goudsmit.) 

Turning to quantum mechanics, relativistic behavior is 
here introduced as a perturbation to the non-relativistic 
Coulomb potential. To discuss this type of precession in 
quantum mechanics, the first requirement is a quantum 
mechanical orbit. Eigenstates of the unperturbed Hamil- 
tonian are obvious candidates, but the standard states 
\n,£,m) are certainly not elliptical. However, for this 
system there exist generalized coherent states which are 
localized along bound classical trajectories of arbitrary 
eccentricity 0. These states are formed of the eigen- 
states pertaining to a single energy level and are hence 
themselves eigenstates of the Hamiltonian, stationary in 
time. With the introduction of a special relativistic per- 
turbation, the states will change in time, and in the limit 
of large quantum number, should precess at the rate pre- 
dicted by special relativity. As shown below, this is in- 
deed the case. 
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The perturbation to be considered is given by 
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so that the first order Raleigh-Schrodinger perturbative 
correction to the energies is given by 
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in which m is the electron mass, c is the speed of light, 
Z is the atomic number, a = e 2 /Tic is the fine structure 
constant, e is the electron charge, n is the total quan- 
tum number, and £ is the quantum number pertaining 
to angular momentum, the eigenvalues of L 2 given by 
h 2 i(£+ 1). The perturbation Eq. fl2.3| ) is obtained by ex- 
panding the special relativistic expression for the kinetic 
energy and retaining the next term after the Newtonian 
term. No spin-orbit coupling is considered in parallel to 
the classical theory to which these calculations will be 
compared. 

The classical peri od may be extracted from the quan- 
tum spectrum Eq. fl2.2| ) as follows. The expansion of the 
energy eigenlevels about n = (n) is given by 
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Now assume that the system is in some state 
with average total quantum number (n) and uncertainty 
(An) 2 = (n 2 ) — (n) . If the quantum time evolution is 
expressed in terms of the eigenstate expansion, 



\m) = J2< 



-iE„t/h 



(2.6) 



II. ESTIMATING THE RATE OF PRECESSION 

In the non-relativistic theory, the unperturbed Hamil- 
tonian is given by 



the first term in Eq. (2.5) leads to an overall phase factor 
which may be dropped, and the second term will lead to 
phase factors that are integer multiples of 27r when 
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This gives the classical period of a classical trajectory 
with energy E^ exactly, a circumstance closely con- 
nected to the success of Bohr-Sommerfeld quantization 
in hydrogenic atoms. Also, for hydrogenic wave func- 
tions, radial expectation values are given by 
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Substitution of this value into T c i yields 
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a quantum rendering of Kepler's third law. 

If the expansion Eq. (^?j) terminated with the linear 
term, then periodic motion would continue indefinitely. 
However, the quadratic and higher terms contribute a 
dephasing influence and the approximation of periodic 
behaviour eventually breaks down. To estimate the size 
of this dephasing influence, evaluate the quadratic term 
at the edge of the distribution in n (n = (n) ± An) at 
t = Tci to yield the "test quantity" 
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Hence, in the special case of a state localized in position 
and momentum, Ehrenfest's equations indicate an initial 
trajectory which follows a classical trajectory. In the 
additional case that 6<p -C 2ir, the state achieves at least a 
full period of behavior approximating classical behavior. 
These considerations are similar to those used employed 
in describing wave function revivals |^| [J . 

Now treating the relativistic perturbation to kinetic 
energy Eq. ( |2.4D in a similar manner^ the expansion about 
£ — {£) of the perturbation Eq. (|]J) is given by 
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Suppose the system is now in a state composed of eigen- 
states pertaining to the same n, degenerate in the unper- 
turbed syste m. Then, the term independent of (£ — (£) ) 
in Eq. ( 2.11 ) leads to overall phase factors, and the linear 
term leads to integer multiples of 2n when 



With the additional assumption that the system is in 
a state built up upon a particular classical orbit (such 
states are described in the next section), the evolution 
which occurs on a time scale T p derived above must be of 
the same nature as classical precession due to relativity, 
that is, in-plane rotation of the state about the origin 
with a period of T p . If this is the case, then the angle of 
rotation per classical period T c i is given by 
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Co ntinu ing the analogous arguments whi ch lead to 
Eq. ( 2.1C ), the quadratic term of Eq. fl2.1l| ) evaluated 
at t = T p at the edges of the distribution in I yields the 
test quantity 
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Again, if 5<j) is small, then the state at t — T p will ap- 
proximately resemble the state at t = 0. This condition 
is best satisfied by Rydberg states with large angular mo- 
mentum. 

From Bergmann JtJ , the classical rate of perihelion pre- 
cession according to special relativity is given byQ 
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measured as a change in angle per classical period of 
the vector pointing from perihelion to aphelion. In this 
expression, G is the gravitational constant, M is the 
(large) mass of the source of the gravitational field, m 
is the (small) mass of body in orbit, and L is the orbital 
angular momentum of the o rbiti ng body. To compare 
the quant um prediction Eq. (2.14) with the classical rate 
Eq. ( 2.16| ), replace the strength of the Coulomb potential 
(Ze 2 ) with that of the Kepler problem (GMm), and note 
that the total angular momentum squared will be close 
to h 2 ((£) + \) 2 given the above assumptions that (£) is 



(^■ll) large and (A£) 2 is small: 
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Thes e su bstitutions render the quantum estimate 

Eq. (2.14) identical to Eq. ( 2.1 6| ) as anticipated. 
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This may be written in terms of T c /, yielding 
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x This rate is never observed in nature since those objects 
whose orbits are observed to precess are invariably gravi- 
tationally bound, requiring a general relativistic treatment. 
Such a treatment yields a rate, observed in the case of Mer- 
cury §, of Slo = 6-KG 2 M 2 m 2 /c 2 L 2 . 



Note that T p > T cl which follows from E^\ < . 
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III. SO (4) COHERENT STATES 

The generalization of coherent states due to Pcrclo- 
mov H,[10| is useful in this context. This generalization 
rests upon the group structure of the system in ques- 
tion. In the case of the hydrogen atom, the dynamical 
group is SO (4, 2) 1 11 - 13 j. For the present purposes, it is 
not necessary to engage the entire group; the degeneracy 
group SO (4) contains sufficient structure. Treatments 
of coherent states of this description are existent in the 
literature [§§]. 

The realization of SO (4) which describes the degen- 
eracy of the hydrogen atom spectrum is given by the 
elements of the angular momentum operator and those 
of the scaled Laplace-Runge-Lenz vector. Given a classi- 
cal orbit, the classical version of the Laplace-Runge-Lenz 
vector is proportional in magnitude to the eccentricity 
and is aligned parallel to the major axis. The corre- 
sponding quantum operators in atomic units are given 
by H 



L = r x p, 



(3.1) 
(3.2) 



All of these operators commute with the Hamiltonian, 
leading to conservation of these quantities under time 
evolution. Conservation of angular momentum is to be 
expected in this (spherically symmetric) potential, but 
conservation of the A (a simple derivation of which is 
given by Wulfman juj) is a unique property of the non- 
relativistic Coulomb potential. The classical interpreta- 
tion of the invariance of the Laplace-Runge-Lenz vector 
is that elliptical orbits do not precess. 

These operators satisfy the commutation relations 

[Lj,Lk] = iejkiLi, [Aj,Ak] = iejkiLi, [L ,A k ] = ie jk iAi. 

(3.3) 

These six operators may be decoupled into two groups of 
three operators via 
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M j = -(L j +A j ), N 3 = -(L 3 -A 3 ), 
which commute according to 



(3.4) 



[M 0l M k ] = ie jH Mi, [Nj , N k ] = ie jk iN h [Mj , N k ] = 0. 

(3.5) 

In terms of these new operators, it is clear that SO(4) = 
SO(3)®SO(3). 

The group SO (4) has two Casimir operators, given by 



Ci = L 2 + A 2 = 2{M 2 + N 2 ), 
C 2 = t-A = M 2 -N 2 . 
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In the hydrogenic realization of this group, Ci — so that 
quantum mechanically and classically, angular momen- 
tum is perpendicular to the Laplace-Runge-Lenz vector. 
The second Casimir operator also indicates that the di- 
mensions of the irreducible representations of the SO (3) 
generated by the M k and by the N k are equal (say, to 
n = 2j + l) so that the dimensions of the relevant unitary 
irreducible representations of SO (4) are n 2 , the famous 
degeneracy of the hydrogen atom energy spectrum. Thus 
the first Casimir operator is equal to 



Ci - 4j(j + 1) 
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representing a constraint on the sum (L 2 + A 2 ). 

Turning to the SO (4) coherent states, Since SO (4) = 
SO(3) <g> SO(3), the SO(4) coherent states may be ex- 
pressed as the direct product of two SO (3) coher- 
ent states which are themselves standard in the litera- 
ture [p|[l0 16|. The SO(3) coherent states are given by 
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parameterized by the complex valued £. In these states, 
expectation values of the angular momentum operators 
are given by 



T ^(Re(C),-Im(C),i(|C| 2 
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with J standing for M or N as the case may be. With 
these expressions in mind, the SO (4) coherent states are 
given by 



hCi,C 2 ) = |i,Ci>li,Ca>. 



(3.11) 



Here, n = 2j + 1, and £i and ^2 parameterize the SO (3) 
coherent states pertaining to M and N respectively. 
The expectation values of the angular momentum and 
Lapl ace-R unge-Lenz operators may be regained through 
Eq. (|3.10|) and the relations L = M + N an d A — M-N. 
(The substitution of Eq. |T|) into Eq. ( [3.1l| ) yields a 
double sum over the direct product states |j, mi)\j, m?). 
These states may be calculated by relating them to the 
standard hydrogenic eigenstates \n, £, m) via the Clebsch- 
Gordon coefficients.) 

For the purpose of visualization, spherical symmetry 
permits setting (L) parallel to the z-axis, achieved by 
setting (2 = — Ci- Since L • A = 0, the vector (A) may be 
oriented parallel to the x-axis accomplished by setting 
the imaginary parts of £1 and (2 equal to zero. These 
identifications reduce the problem to the variation of a 
single real parameter, say 77 = Re(fi)- In terms of this 
parameter, expectation values are given by 
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With the magnitude of the classical Laplace- Runge-Lenz 
vector being proportional to the eccentricity e of the or- 
bit, the quantum calculation leads to 
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Also in terms of rj, the total angular momentum may be 
expressed as 



(L 2 ) = 2j(j + 1) + 2j 



, 2 rf — 6rj 2 + 1 



(l + V 2 ) 2 
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This expression fol- 

lows from {P) = ((M + N) • (M + N)), Eq. pi), and 
(M • N) = (M) -(N), the last of which since [Mj,N k ] = 0. 




FIG. 1. Eccentricity e versus rj, the parameter to the SO(4) 
coherent state. 



By virtue of their construction, the SO (4) coherent 
states are minimum uncertainty states of rotated 4- 
dimensional angular momentum operators. One there- 
fore expects that for moderate 4-dimensional rotations, 
the states will remain localized in angular momentum 
and in the Laplace-Runge-Lenz vector, that is, localized 
about a particular classical orbit. Hence, given a valu e 
of rj, the approximate geometry is gi ven by Eq. (3.13) 
with a semi-major axis given by Eq. 



IV. PRECESSING COHERENT STATES 



The value of the test quantity 8(f) given by Eq. ( 2.15| ) 
is expressible in terms of the SO (4) coherent state pa- 
rameter rj. With expectation values taken in an SO (3) 
coherent state, 



(AJ 3 ) 2 = 2j 
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and 
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commute. With these in mind, and 
from the denominator of 8<f> since (£) is 



= 2tt- 



1 



(4.3) 



Therefore, in the sense that 8(f> = 0, the precession is 
best for 77 = or as 77 — » 00. In both cases, though, from 
Eq. (3.13), the resultant orbits will be circular and no 
precession can be observed: There is a certain tradeoff 
between the observability and "coherence" of the preces- 
sion measured by S(f>. In f act, i t is convenient to express 
8(f> in terms of e from Eq. ( 3.1 3| ) so as to connect this test 
with a more physical or geometrical quantity: 



2ne 2 (2- e 2 ±2V T^ 
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where the positive sign is appropriate for rj > 1 and the 
negative for 77 < 1. As shown in Figure |^, as the eccen- 
tricity increases, the degree to which the wave functions 
remain assembled decreases. 




FIG. 2. Test parameter 8<j> versus eccentricity e. 

For some examples of precessing coherent states, exam- 
ine Figure |[ This figure depicts a mesh plot of an SO (4) 
coherent state at t — 0, and a sequence of overlayed 
contour plots showing states precessed from t = to 
t = jT p . To appreciate the physical scale of these simu- 
lations, the field of view in all cases is 4.23 /im across and 
with states in the 141st energy level, the classical period 
is 4.25 x 1CT 10 seconds. In Figure |(b), with e = 0.385, 
the precession period is 0.266 seconds or 6.26 x 10 8 clas- 
sical periods, which on atomic scales, is a very long time. 
In the subsequent images, the rate of precession is larger 
owing to the larger eccentricity so that the precession pe- 
riod of the state depicted in Figure ||(d) is 0.164 seconds. 
This interval is still much longer than the longest times 
Rydberg states are observed in experimental setups, typ- 
ically 1 nanosecond ||. 
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a more involved classical analysis. In particular, a spe- 
cial relativistic perturbation to non-relativistic quantum 
mechanics leads to an agreement with classical special 
relativity in a large quantum number limit. A similar 
notion of agreement may be found in the work of McRae 
and Vrscay JlTj who have studied correspondence be- 
tween quantum and classical perturbation schemes. In 
their work, as well as in the present paper, it transpires 
that the simplest route to the determination of a classi- 
cal perturbation may be through the classical limit of a 
quantum perturbation. 




(d) 



v 



(if 


fey 







-4 -2 2 4 -4 -2 2 4 

X X 

FIG. 3. Some examples of precessing SO(4) coherent 
states: (a) In the 141st energy level, at t = 0, with r\ = 0.2 
(or e = 0.385) on the x-y plane with amplitude proportional 
to | (r 1 141, £i, £ 2 ) | 2 - (b) The same state as (a), with posi- 
tion plotted in 10 4 atomic units, and the origin located at the 
+ symbol. The state at t — is oriented horizontally, and 
t = jT p oriented vertically. The solid lines depict a single 
contour on the quantum wave function, and the dashed lines 
depict the classical orbit precessed according to the special 
relativistic prediction, (c) The same as (b), but with r\ = 0.3 
(or e = 0.550). (d) The same as (b), but with 77 = 0.4 (or 
e = 0.690). 

As apparent from the images, the quantum rate of pre- 
cession agrees with the classical special relativistic pre- 
diction. In the case of Figure ||(b), though only shown 
until t = jT p , the state remains localized on the ellipse 
up to the full period of the precession. Not surprisingly, 
for increased eccentricities, the degree of "coherence" of 
the state decreases, leading to the decay of the state as 
depicted in Figure [|(d) after only a quarter precession 
period. 



V. CONCLUSIONS 



The perturbation Eq. (2^3) was chosen so that a di- 
rect comparison could be made between the quantum 
mechanical calculation and the special relativistic calcu- 
lation due to Sommerfeld. Therefore, spin effects were 
ignored which are equal in size to the kinetic energy 
perturbation. This means that the quantities calculated 
here are not predictions to be tested in the laboratory. 
The purpose has been to show that a simple analysis 
of a quantum perturbation can reproduce the results of 



VI. ACKNOWLEDGMENTS 

This work was supported by the Natural Sciences and 
Engineering Research Council of Canada, and the Na- 
tional Science Foundation. The author would also like to 
acknowledge useful discussions with Prof. N.W. Ashcroft 
of the Department of Physics, Cornell University, and 
Profs. E.R. Vrscay and J. Paldus of the Department of 
Applied Mathematics, University of Waterloo. 



[1 
[2 

[3; 
[5 

ir 

[§ 

[9 
[10 

[11 



A. Sommerfeld. Atombau und Spektrallinien. F. Vieweg 

und Sohn, Braunschweig, Germany, 1921. 

J.C. Gay, D. Delande, and A. Bommier. Atomic quantum 

states with maximum localization of classical elliptical 

orbits. Phys. Rev. A, 39(12):6587-6590, 1989. 

I.S. Averbukh and N.F. Perclman. Fractional revivals: 

Universality in the long-term evolution of quantum wave 

packets beyond the correspondence principle dynamics. 

Phys. Lett. A, 139(9) :449-453, 1989. 

M. Nauenberg. Autocorrelation function and quantum 

recurrence of wavepackets. J. Phys. B, 23:L385-L390, 

1990. 

J. A. Yeazell, M. Mallalieu, and Jr. C.R. Stroud. Obser- 
vation of the collapse and revival of a Rydberg electronic 
wave packet. Phys. Rev. Lett, 64(17):2007-2010, 1990. 
M.G.A. Crawford. Temporally stable coherent states in 
energy degenerate systems: The hydrogen atom. Phys. 
Rev. A, 62(l):012104-l-7, 2000. 

G. Bergmann. Introduction to the Theory of Relativity. 
Prentice-Hall, Englewood Cliffs, NJ, 1942. 
CM. Will. Theory and Experiment in Gravitational 
Phystcs. Cambridge University Press, Cambridge, revised 
edition, 1993. 

A.M. Perelomov. Coherent states for arbitrary Lie 
groups. Commun. Math. Phys., 26:222-236, 1972. 
A.M. Perelomov. Generalized Coherent States and Their 
Applications. Springer- Verlag, London, 1986. 
A.O. Barut and G.L. Bornzin. SO(4,2)-Formulation of 
the symmetry braking in relativistic Kepler problems 
with or without magnetic charges. J. Math. Phys., 
12(5):841-846, 1971. 



5 



[12] J. Paldus. Dynamical groups. In G.W.F. Drake, editor, 
Atomic, Molecular and Optical Physics Handbook. AIP 
Press, Woodbury, New York, 1996. 

[13] B.G. Adams, J. Cizek, and J. Paldus. Lie algebraic meth- 
ods and their applications to simple quantum systems. 
In Advances in Quantum Chemistry, Vol 19. Academic 
Press, New York, 1988. 

[14] C.E. Wulfman. Dynamical groups in atomic and molec- 
ular physics. In E.M. Loebl, editor, Group Theory and 
Its Applications, volume II. Academic Press, New York, 
1971. 

[15] J.R. Klauder and B.-S. Skagerstam, editors. Coher- 
ent States: Applications in Physics and Mathematical 
Physics. World Scientific, Singapore, 1985. 

[16] W.-M. Zhang, D.H. Feng, and R. Gilmore. Coherent 
states: Theory and some applications. Rev. Mod. Phys., 
62(4):867-927, 1990. 

[17] S.M. McRae and E.R. Vrscay. Perturbation theory and 
the classical limit of quantum mechanics. J. Math. Phys., 
38(6):2899, 1997. 



G 



